Pinning and transport of cyclotron (Landau) orbits by electromagnetic vortices 
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Electromagnetic waves with phase defects in the form of vortex lines combined with a constant 
magnetic field are shown to pin down cyclotron orbits (Landau orbits in the quantum mechanical 
setting) of charged particles at the location of the vortex. This effect manifests itself in classical 
theory as a trapping of trajectories and in quantum theory as a Gaussian shape of the localized wave 
functions. Analytic solutions of the Lorentz equation in the classical case and of the Schrodinger 
or Dirac equations in the quantum case are exhibited that give precise criteria for the localization 
of the orbits. There is a range of parameters where the localization is destroyed by the parametric 
resonance. Pinning of orbits allows for their controlled positioning. They can be transported by the 
motion of the vortex lines. 
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INTRODUCTION 



The transverse motion of charged particles in a con- 
stant magnetic field is fully delocalized. The classical 
trajectories of particles projected on the plane perpen- 
dicular to the field are circles which can be moved around 
without changing the particle energy. The wave functions 
of a particle in quantum mechanics exhibit the same be- 
havior. If vb(r) is a solution of the Schrodinger or the 
Dirac equation with the energy E, then a shifted wave 
function exp(ier-(?"o x B)/2h)ip(i — ro) is also a solution 
with the same energy. In the present paper we shall show 
that electromagnetic beams with vortex lines localize the 
classical and quantum states at the position of the vor- 
tex and when the vortex line is moved the orbits will fol- 
low. In order to describe these properties quantitatively 
we shall use new exact solutions of the equations of mo- 
tion obtained in the presence of an electromagnetic field 
comprising a constant magnetic field and a wave with a 
vortex line. This paper extends our earlier analysis by 
adding a constant magnetic field. This extension is signif- 
icant because the addition of a constant field introduces 
the third parameter characterizing the electromagnetic 
field. From the three parameters (the wave frequency 
and two field amplitudes) we can construct two dimen- 
sionless parameters and the space of distinct solutions, 
as compared to those given in 0, becomes much richer. 
An extensive analysis of the exact solutions of relativis- 
tic equations of motion can be found in a monograph by 
Bagrov and Gitman 2]. However, they overlooked the ex- 
istence of the solutions described in this paper. A related 
problem — the motion of a particle in a constant mag- 
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netic field and a plane electromagnetic wave — has been 
treated in detail by Roberts and Buchsbaum Q a long 
time ago. The crucial difference between these two cases 
is that the translational symmetry in the plane perpen- 
dicular to the magnetic field is broken by the presence 
of an electromagnetic vortex and the problem becomes 
truly three-dimensional. 

The aim of the present paper is not only to present new 
analytic solutions of the classical and quantum equations 
of motion in some realistic configurations of the electro- 
magnetic field but also to describe a universal confining 
mechanism of charged particles that might have exper- 
imental applications. This mechanism acts for all elec- 
tromagnetic waves with a definite angular momentum. 
Such waves are characterized by the presence of an elec- 
tromagnetic vortex. At the vortex the field vanishes and 
near the vortex the electric and magnetic field vectors ex- 
hibit a characteristic behavior that leads to the trapping 
of particles. The trapping of atoms by electromagnetic 
vortices is a well established phenomenon It em- 

ploys the dipole force pushing atoms in the direction of 
a decreasing field amplitude in a laser beam that is blue- 
detuned from the relevant atomic transition. The mech- 
anism for trapping charged particles by electromagnetic 
vortices is completely different. It employs the rotation 
of the electric and magnetic field vectors near the vortex 
line. 

This paper contains a complete, classical, and quantum 
mechanical description of the motion of charged particles 
in a combination of an electromagnetic wave with an em- 
bedded vortex and a constant magnetic field. Our treat- 
ment is fully relativistic but we also present the results in 
the nonrelativistic limit. In Sec. II we introduce a model 
electromagnetic wave that approximates very well, near 
the vortex line, all solutions of Maxwell equations with 
the lowest angular momentum. In Sec. Ill we give gen- 
eral solutions of the Lorentz equations in such a wave 
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and exhibit a mechanism that confines the particle in 
the vicinity of the vortex line. In Sec. IV we show how 
trapped particles can be transported by moving vortex 
lines. In Sec. V we carry our analysis to the quantum 
theory presenting exact solutions of the Dirac equation 
in the presence of a wave with a vortex and the magnetic 
field. Finally, we present solutions of the Schrodinger 
equation describing the transport of Landau orbits. 



II. THE SIMPLEST ELECTROMAGNETIC 
BEAM WITH A VORTEX LINE 



The notion of a vortex line of the electromagnetic field 
can be traced back to an old paper by Rainich A 
detailed analysis with many examples has been recently 
given by us in Ref. 9]. In thia paper we shall use only 
a very special case of the vortex line — the one that 
is found in null solutions of the Maxwell equations [To| . 
In reality, vortex lines of the electromagnetic field are 
most commonly found in beams carrying angular mo- 
mentum [ill Ha . [l3j . Such beams can now be produced 
with the use of several methods: computer generated 
holograms, axicons, spatial light modulators, and biaxial 
crystals, etc.. The best known examples of such beams 
are the Bessel beams and the exact Laguerre-Gaussian 
(LG) beams. These are solutions of the Maxwell equa- 
tions obtained by separating the variables in cylindrical 
coordinates 0] . A convenient description of these beams 
in terms of just one complex function is obtained with the 
use of a complex vector 



F(r, t) = E(r, t) + icB(r, t), 



(1) 



which we named the Riemann-Silberstein (RS) vector 
[THl Hfjj . The RS vector obeys the complex form of 
Maxwell equations 



d t F(r,t) = -icV x F(r,t), V-F(r,t) = 0. 



(2) 



We shall use a representation of F in terms of just one 
complex function 



F x = (d x d z + -dyd t )x, 

1 

F y = (d y d z - -d x d t )x, 



-{dl + dl) X . 



(3a) 

(3b) 
(3c) 



where x is an arbitrary complex solution of the wave 
equation 



5 #-A)x(r,t) = 0. 



(4) 



An equivalent representation of the electromagnetic field 
expressed in terms of two real solutions of the wave equa- 
tion, instead of a single complex solution of the wave 
equation, has been given by Whittaker [HI- The func- 
tions x( r jt) f° r the Bessel beams and for the Laguerre- 
Gaussian beams are most easily expressed (cf. |18|1 in the 



cylindrical coordinates (p, cf), z). The Bessel beams are la- 
beled by the following parameters ("quantum numbers"): 



the transverse wave vector k± = yj fc| + ky, the (dimen- 

sionless) angular momentum m, the wave vector k z in 
the direction of propagation, and the hclicity a = ±1 



XB(p,(j),z,t) = e 



— ia{ijj}.t— k z z— m> 



«J m (fcxp). (5) 



The Laguerre-Gaussian beams are defined as 



XLG(p,<t>,z,t) 



a(t+) 



n+m+l 



■ exp 



a(t+) 



(0) 



where t± = t ± z/c and a(t + ) = I 2 + iac 2 t + /uj, and L™ 
is the Laguerre polynomial. The natural number n gives 
the number of zeros of the polynomial but otherwise has 
no direct physical meaning. The parameter I determines 
the transverse size (waist) of the LG beam. In the limit, 
when k± — ► for Bessel beams and I — * oo for LG beams, 
both functions (JSJ) and © reduce (for m > and apart 
from numerical prefactors) to the same solution of the 
wave equation 



(7) 



This choice of x leads to the following RS vector 

F(r, t)=a{k + iy) (x + iy) m ^ e~ lau)(t - z ' /c \ (8) 

where we again disregarded some numerical prefactors. 
The RS vector (jSJ) is a good approximation to the 
Laguerre-Gaussian and Bessel beams near the z axis, 
when p <C I or pk± <C 1 . It belongs to a family of the so- 
lutions of the Maxwell equations describing vortex lines 
riding atop a null electromagnetic field (lffi 

The simplest solution, the one with a unit vortex 
strength, is obtained by choosing m = 2 in Eq. JSJ. In 
this case, the electric and magnetic field vectors are 

E(r,t)=aBLj(f(r,t),g(r,t),0), (9a) 

B(r, t) = a— (-g(r, t), /(r, t), 0) , (9b) 
c 

where we introduced the field amplitude expressed in 
terms of the magnetic field strength and 

f{r,t) — x cos ui(t — z/c) + aysinui(t — z/c), (10a) 
g(r, t) — ax smu(t — z/c) — ycosuj(t ~ z/c). (10b) 

This simple, linear dependence of the field vectors on 
x and y will enable us to find explicit solutions of the 
equations of motion. The expressions © and (|10|l show 
that the sign factor a determines the helicity — the sense 
of rotation of the pair of vectors E and B. At each point 
this pair of vectors rotates with the wave frequency u> 
but, in contrast to the circularly polarized plane wave, 
the orientation is not the same in the whole plane but 
the pair rotates as we move around the vortex line (see 



3 




FIG. 1: This figure shows the orientation and also the 
strength of the electric and magnetic field vectors at t = 
and in the z — plane for the solution of the Maxwell equa- 
tions with a vortex line described by Eq.©. The point in the 
middle represents the position of the vortex line. All plots in 
this paper were generated with the help of MATHEMATICA 





Fig-C]). I 11 our case, since the vortex has a unit strength, 
the full rotation angle is 2tt. 

We have shown recently [| that there exist analytic 
solutions in the classical and quantum theory describing 
the motion of charged particles in the presence of the 
electromagnetic beam [Eq|5]. Exact solutions exist be- 
cause the transverse motion (in the plane perpendicular 
to the direction of the beam propagation) almost sepa- 
rates from the longitudinal motion. The separation is not 
complete but the longitudinal motion influences only the 
values of the integration constant appearing in the equa- 
tions of motion describing the transverse dynamics. The 
exact solvability is not destroyed by the presence of an 
additional constant magnetic field aligned with the beam 
direction. The presence of this field adds a new dimen- 
sion to the space of control parameters and significantly 
enriches the dynamical behavior of the system. 

The electromagnetic field © may seem unrealistic 
since its amplitude increases without bound with the 
growing distance from the z axis. However, we shall show 
that the particle is often confined to a region close to the 
z axis and is subject practically the same field as that of 
much more realistic Laguerre-Gaussian or Bessel beams. 
This expectation has been recently fully confirmed in 0] 
where we have shown that numerical solutions of the 
Lorentz equation in the presence of Bessel beams near 
the vortex line exhibit the same behavior as the analytic 
solutions found in [| for the model wave ©. 



III. CLASSICAL MECHANICS 

We shall give first a fully relativistic description of the 
motion because the nonrelativistic approximation can be 
easily obtained from the relativistic solution by taking 
the limit c — > oo. 



A. Relativistic motion 

The Lorentz equations of motion 

me(T)=en(r)}Ur), (11) 

that govern the motion of a relativistic particle moving 
in the presence of the electromagnetic wave [Eq.©] and 
a constant magnetic field Bq — (0,0, Bo), aligned with 
the beam direction, can be written in the form 

£ = ouj c uj /(£, n, C, 6) (§ - C/c) + w fi, (12a) 
ri = au c ujg(t,r,,(,0)(6-(/c) -wqI (12b) 
C = <^ [Zf(Z, V, 0) + r,g(Z, V, 0)}, (12c) 

r J = a^tf(Z, V , C, 9) + r,g(£, r,, 9)}, (12d) 
c 

where £(r), r](T), C(t) are the space components and 9(t) 
is the time component of the four- vector £ M (r) . To save 
space we dropped the dependence on r. The dots de- 
note derivatives with respect to the proper time r. The 
two cyclotron frequencies lo c = eB/m and ujq = eBo/m 
measure the intensity of the wave and the strength of the 
constant magnetic field, respectively. They can be either 
positive or negative. The negative value of loq means the 
reversal of the magnetic field direction. 

From now on, we shall assume that a = 1 because 
the equations of motion (|12fl with the functions / and 
g defined by Eq. i|10[l are invariant under the following 
transformation 

a -> -cr, £ < — ► rj, uj -> -w . (13) 

Therefore, we do not obtain anything essentially new by 
considering both signs of a, since to every trajectory for 
a = — 1 there corresponds a trajectory for a = 1 in the 
reversed constant magnetic field, differing only by a re- 
flection with respect to the x — y plane. 

The Lorentz equations (|12f> arc nonlinear but they can 
be effectively linearized owing to conservation laws. In 
doing so, we shall follow the procedure developed in 
Ref. 0| . The first constant of motion 8 is obtained 
by subtracting Eq. i|12c[l from Eq. (|12d(l which gives 
9 — Ql 'c = or 9 — Q) 'c = consti. Thus, £ is a conserved 
quantity 



£ = 9 - C/c = yfl + (Z 2 + V 2 + C 2 )/c 2 - C/c (14) 
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Apart from the factor mc 2 , the constant £ is the light - 
front energy — the conjugate variable to 4+ = t + z/c 



£ 



1 — Vz/c \J m 2 c 4 + p 2 c 2 



PzC 



^Jl-v 2 /c 2 



mc* 



(15) 



We shall assume that we start counting the proper time 
t in such a way that 9(0) — C(0)/ c = 0. Under this 
assumption, Eq. l |T3j l integrated with respect to r yields 



C/c = £> 



(16) 



Hence, the proper time is simply proportional to the 
light-front variable and this will enable us to separate 
the transverse motion from the longitudinal motion. 

The second constant of motion is obtained by subtract- 
ing from Eq. (|12cl) the sum of Eq. I|12afl multiplied by £/£ 
and Eq. (|12bfl multiplied by f\j£. Since the right hand 
sides cancel, we obtain 



c-— (i 2 

^ 2c£ ^ 



= const 2 = ^(i-£). (17) 



The last expression is obtained by rearranging the for- 
mula 1 + (£ 2 + fj 2 + C 2 )/c 2 = (£ + C/ c ) 2 obtained from 
Eq. O. 

Owing to Eq. l|14f) . the longitudinal motion affects the 
transverse motion only through the value of £ and since 
it is a constant of motion the equations for £ and r\ may 
be solved first. The integration of Eqs. (|12afl and l|12b|) 
is made simpler by introducing a complex combination 
of the coordinates Z = £ + irj. The equation of motion 
for Z obtained from Eqs. |Q reads 



Z — uj c QZ*e 



•fir 



iuj Z, 



(18) 



where in addition to the two cyclotron frequencies ui c 
and wo associated with the strength of the wave and the 
constant magnetic field we defined the third frequency 
ft — an effective wave frequency as seen by the moving 
particle. Thus, the dynamics of the particle is governed 
by the three frequencies 



il = LOG, LU C = , LOq = — • 

m m 



(19) 



The form Ze~ mT l 2 = io c n(Z e - inT / 2 )* - iuj Ze- inT ' 2 of 
Eq. (|T%)l suggests a transformation to the frame rotating 
(in proper time) with the angular velocity ft/2 around 
the z-axis. This amounts to replacing £ and 77 by the 
new variables 



a(r) = f (t) cos(fir/2) + r/(r) sin(fir/2), 
Pi T ) = -£(t) sin(fir/2) + t)(t) cos(fir/2), 

or in complex notation 



(20a) 
(20b) 



(21) 



Upon substituting this expression into Eq. (|18fl . we get 
rid of an explicit r dependence and obtain the following 
equations with constant coefficients for a and j3 

a = (fi 2 /4 + Qujq/2 + Quj c )a + (fl + lj q )/3, (22a) 
= [n 2 /A + fiwo/2 - fkj c )/3 - (ft + w )d. (22b) 

The mutual relationships between the three frequencies 
ft, io c, and loo determine whether the transverse motion 
is localized near the z axis or is unbounded. 

A full description is most easily done in the Hamilto- 
nian formalism. The equations of motion (|22|l can be 
obtained from the following Hamiltonian 

„ Pl+Pl , (a - b)a 2 + (a + b)(3 2 
H = — h m 



2m 2 
-wtapp - (3p a ), 



where 



fllOr 



ft 



- w 



(23) 



(24) 



This Hamiltonian (for luq = 0) differs by a canonical 
transformation from the one considered in Ref. 0. We 
found this new form of the Hamiltonian more convenient 
than the previous one. The canonical equations of mo- 
tion can be written in the following matrix form 



where 



and 



i—X = X-M, 
dt 



X = (a, (3, p a , pn ) 



M 



( 

iw 

i/m 




—iw 



i/m 



-im(a - 





b) 





-im(a 4 
—iw 




b) 



(25) 



(26) 



(27) 



We have absorbed i into the definition of the matrix M 
to make its eigenvalues equal to the (real) characteris- 
tic frequencies. These are the same equations of motion 
(only with different values of the coefficients) that govern 
the dynamics of particles in the mechanical model of the 
Paul trap [2(], |2y, Trojan asteroids in the Sun- Jupiter 
system, and Trojan wave packets of Rydberg electrons 
in an electromagnetic wave \22\ or in a molecule with an 
electric dipole moment [23|. As a matter of fact, every 
quadratic Hamiltonian, by an appropriate linear canoni- 
cal transformation, can be reduced to the form l|23|) . 
The characteristic frequencies in the present case are 



ft j 



ftr± 



r± = rV(l + M) 2 + V? ± 4^ 2 + M 2 (l + m) 2 /4 , (28) 
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where we used f2 as a yardstick to measure all frequencies, 
i.e. 



p, = u>o/Q, v = u; c /f2. 



(29) 



The description of the motion with the use of two dimen- 
sionless parameters /i and v is very convenient because 
it exhibits a self-similarity inherent in our problem. The 
dimensionless ratios fi and v provide a unified description 
of different physical situations. For example, a particle in 
a 2.45 GHz microwave oven will show the same behavior 
as a particle in an optical wave of a 505 nm blue laser 
provided we increase the amplitude of the wave and the 
strength of the magnetic field by a factor of 2.4210 5 . An 
increase or decrease of all three frequencies by the same 
factor does not change the character of the motion but 
it results in a decrease or increase of characteristic space 
and time scales. We shall take this fact into account and 
express the distances in units of the angular wavelength 
X = c/uj and the time in units of an angular period 1/ui. 
These units of length and time are used in all the figures 
in this paper. In other words, we shall use the units in 
which lu — 1 and c = 1 . 

Bounded oscillations occur only when both character- 
istic frequencies f2± are real and this means that 



\»\<\ 



1 

2+^ 



(30) 



The mode amplitudes a± and a* ± (classical counter- 
parts of the annihilation and creation operators) are 
found by solving the eigenvalue problem for the matrix 
M appearing in Eq. I|25|) 



N 



1/2 



P/3 



—let 
1 



N 



1/2 



Pa 

t— [ PS3 



+iet-+ p a 



mQs- 
1 + M ' 

' r+/T / 

mf2s+ 

( 

1 + M 

-P 



1 



(31a) 



(31b) 



where to save space we introduced the following functions 
of jjL and v 



v/^ 2 +M 2 (l + ^) 2 /4±^ 
v/(l + /i) 2 + 2s±, 



2s. 



t_± = VI(i + m) 2 

e = sgn(l + /i), 
and the normalization factors are 

N± = 4mfi±(s+ + s_). 



(32a) 

(32b) 

(32c) 
(32d) 

(33) 



This normalization guarantees that a± and a± have the 
canonical Poisson brackets 



{a±,a±} 



(34) 



The formulas l|31|l are not valid when fx = — 1. However, 
in this special case the equations of motion (12211 describe 
just two uncoupled harmonic oscillators and the deter- 
mination of eigenmodes is very simple. The Hamiltonian 
[Eg. 123(1 ] expressed through the amplitudes a± reads 



H = VL, 



sgn(l/2 + fijfl- a*_a-. 



(35) 



The minus sign in the diagonal form of the Hamiltonian, 
indicates that for 1/2 + (i > the oscillations in the 
transverse plane are bounded due to the Coriolis force — 
a characteristic feature of the Paul trap or Trojan aster- 
oids and Trojan electrons. We would like to stress that 
despite its quadratic form the Hamiltonian of the trans- 
verse motion does not correspond to a simple harmonic 
oscillator in two dimensions because the frequencies of 
oscillations fi± are not fixed once and for all but they 
depend on the initial velocities through the parameter £ . 
The regions of stable motion in the pi,v plane are shown 
in Fig. 



FIG. 2: Regions of stability in the [w plane are shown as 
two shaded areas. In the area on the left hand side we have 
a normal harmonic oscillator while in the area on the right 
hand side we have the Trojan regime — stable oscillations but 
the Hamiltonian is not positive definite. 

The general solution for £(r) and j](t) in the labora- 
tory frame is obtained by solving Eq. 1221) in terms of 
eigenmodes and then undoing the rotation in Ea. (l20ll . 
The final expression for the motion of particles in the 
plane perpendicular to the vortex line can be compactly 
written in the same complex form as in 

£(t) + iv(r) = e inr/2 [{iDn + + A) sin(ft + r) 
-(Bk- + iC) sin(f2_T) + (iAn + - D) cos(fi + r) 

+ (C7k_ - iB) cos(rj_r)] , (36) 

but the meaning of A, B, C, -D, and n± is now different, 
namely 



K± 



+ 1/4 + n/2±v 
(1 + n)r± ' 



(37) 



and the constants A,B,C, and D are the following func- 
tions of the initial values of the transverse positions and 
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velocities Having found the solution in the transverse plane, we 

can now just integrate Eq. l|17fl to obtain the following 
_ K y l 1 ~ K - r -) r lQ + So/^ (38a) formula for the longitudinal motion 



(1/2 -K 


- r -) 7 7o 


+ 6>A* 




- r_K_ 


) 


(«+/2 - 






r+ 


— r_K_ 




(1/2 — K 


+?" + )£o 


- j)o/n 


r_ - 


- r + K_ 




(K-/2- 







c= (i/2-« + r +) fo-wn > (38c) 



c(r) = y ( f / dr 2^ (c +?) ) + c(0) - (39) 



£) = ^ ' L -)<*o 10 ' (38d) The integration is cumbersome but elementary and it 

r--r+K-n + leads to 



i (A 2 + d 2 ) ((n 2 + in 2 + )(i + k\) - 8QQ+K+) + (b 2 + c 2 ) ((n 2 + 4n 2 _)(i + kL) - 8nn_«_) 



E- £ 



8c 2 £ 



+ 



+ 



(1 - K 2 + )(fl 2 - Afl 2 + ) [{D 2 - A 2 ) sin(2fi+r) - 2AD(l - cos(2fi+r))] 



32c££L 



(1 - K 2 _)(fl 2 - Afl 2 _) [(B 2 ~ C 2 ) sin(2fi_r) + 2BC(1 - cos(20_t))] 



(« + - K_)(f7 2 - 40 + f7_) + 2(k+k_ - i)nn r 



8c£f2 p 

(k+ + k-)(Q 2 + 4f2 + rj_) - 2(k + k_ + 

8c£ Clm 



(CD - AB) sin(ft p r) - [AC + BD)(l - cos(O p r)) 



(CD + AB) sin(fi m r) - (AC ~ BD)(\ - cos(O m r)) 
I 



C(0), (40) 



where Q p = fi+ + f2_ and Q m = Q + — Sl_. Thus, the 
motion in the z direction is a composition of a uniform 
motion in the proper time r (not in the laboratory time t) 
and oscillations with four frequencies 2Q± and S7 + ±f2_. 
By a fine tuning of the initial conditions we may cancel 
the uniform motion completely and leave only the oscil- 
lations, but any departure from these special values will 
cause a uniform drift. The appearance of two additional 
frequencies f2+ ±f2_ is a new feature that is not found in 
the absence of the constant magnetic field. Note that the 
motion in the transverse plane is not completely decou- 
pled from the longitudinal motion because the effective 
frequency £1, appearing abundantly in the formula (|36|) . 
depends (through £) on the initial velocity in the longi- 
tudinal direction. 

The presence of stability and instability windows is a 
manifestation of the parametric resonance. In the region 
of instability the amplitude of oscillations grows expo- 
nentially. This is essentially different as compared to the 
motion of a charged particle in the circularly polarized 
plane wave and a constant magnetic field. In that case 
there appears a standard resonance, i.e. the amplitude 
grows linearly when the frequency of the driving force 
(the wave) exactly matches the characteristic frequency 
of the system (cyclotron frequency). When the circu- 
larly polarized wave is replaced by a wave with a vortex, 
the set of linear inhomogeneous differential equations be- 



comes a set of homogeneous equations with periodically 
varying coefficients. The driving force disappears but the 
coefficients of the equation become now time dependent 
and there appears a parametric resonance. Its character- 
istic features are: The appearance of the whole regions 
of instability (and not only discrete values as in the case 
of an ordinary resonance), which shrink to just one point 
(u>o = —lu/2) when ui c — > and an exponential growth of 
the amplitude in all regions of instability. 

Three typical stable three-dimensional relativistic tra- 
jectories differing by the choice of initial positions are 
shown in Fig. [31 The trajectories of a particle projected 
on the xy plane for different values of the parameters /x 
and v are shown in Figs. The plots demonstrate 

that the characteristic shapes of the trajectories are very 
sensitive to the changes of these parameters. 



B. Nonrelativistic limit 

The solutions of the equations of motion in the nonrel- 
ativistic regime can be obtained from the relativistic ones 
just by formally taking the limit c — * oo. In this limit, 
the difference between the proper time and the labora- 
tory time disappears. Also, the wave frequency lo and 
the effective frequency fl become equal, because 8 = 1. 
The motion in the longitudinal direction completely de- 
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FIG. 3: Three trajectories obtained for /i = 0.075, v = 0.1, 
and v/c — (0.001, 0, 0), differing only in the initial positions. 

couples from the oscillations in the transverse plane. In 
the nonrelativistic limit c(l/£ — £)/2 — > v z . In this limit, 
the motion in the z direction becomes a free motion de- 
scribed by the formula z(t) = zq + v z t. In contrast to 
the relativistic case, the velocity in the z direction does 
not show any oscillations, v z — v z q. We arrive at the 
same conclusions by solving the nonrelativistic equations 
of motion (a = 1) 



——- = LO c LO (£ COS Lot + 7] sin Lot) + loq — , 

dt z at 



cPr] 

—77 = lo c lo (f sin Lot — -q cos Lot) — loq — , 
dt z 

dt 2 



o. 



(41a) 
(41b) 
(41c) 



which are obtained from the Lorentz equations (|12|l in 
the limit c — > oo. Since the decoupling of the longitudi- 
nal and the transverse motion in the nonrelativistic limit 
is complete, we may easily obtain in this case also the 
solutions for a two-dimensional gas of mutually nonin- 
teracting charged particles kept near the surface by an 
additional confining potential V(z). 



IV. TRANSPORT OF ORBITS BY MOVING 
VORTICES 



Orbits trapped by vortices can be moved around across 
the magnetic field. The simplest case is a uniform mo- 
tion. The electromagnetic field with a moving vortex line 
can be easily obtained from the field with a stationary 
vortex line [Eq.@] by a Lorentz transformation. Assum- 
ing that the motion is along the x axis with velocity v, 
we obtain the following transformed electric and mag- 
netic field vectors 



E v 
~Buj 



cB v 



j(x — vt) cos ip + y sin ip 
j 2 (x — vt) sin (p — 72/ cos ip 
—v{^ 2 {x — vt) cos ip + 7?/ sin Lp)/c 

—"f{x — vt) sin Lp + y cos ip 
7 2 (x — vt ) cos </7 + 72/ sm <P 
v (j 2 (x — vt) sin ip — 72/ cos ip) jc 



(42a) 



(42b) 



where 7 = 1/ ^/l — v 2 /c 2 and (p — to[j(t — vx/c 2 ) — z/c]. 
In Fig. [7] we show the trajectory of a particle that is be- 
ing pulled across the constant magnetic field by the uni- 
formly moving vortex line. Note that this is not the same 
as transforming the whole problem to a moving frame by 
a Lorentz transformation. Such a transformation would 
result also in a change of the constant magnetic field into 
a crossed magnetic and electric fields. 

The transport with a uniform average velocity is just 
the simplest example, but there is a plethora of more 
intricate cases, all of them exhibiting the same behav- 
ior. The simplest method to construct a solution of the 
Maxwell equations with a vortex line undergoing com- 
plicated motions is to superimpose our basic solution 
[Eq.Q] with circularly polarized plane waves moving in 
the same direction, slightly detuned from the frequency 
lo. We shall consider here two cases: a superposition with 
just one plane wave and a superposition with two waves. 
In the first case, the vortex position at a fixed value of z 
will move on a circle and in the second case it will follow 
a trefoil figure. 

In the circular case, we shall choose the Riemann- 
Silberstein vector in the form 

F c (r,t) = Bto (x-My) 

x \(x + iy)e-^ {t ~ z/c] - (xi + ly^e- 1 ^-*'^} , (43) 

where x\ and y\ are two parameters that fix the radius of 
the circle and the initial position of the vortex. The van- 
ishing of the RS vector determines the position (x v ,y v ) 
of the vortex line as a function of z and t 

x v (z, t) + iy v (z, t) = ( Xl + iy^-^-'-^l (44) 



The vortex forms a screw line of the radius \fx\~- ~Yy\ and 
the step 2ttc/(lo— lo\) rotating with the angular frequency 
lo — lo\ around the z axis. The formulas for the electric 
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FIG. 4: The trajectories of a particle projected on the xy plane are shown for different values of the parameter /j, (i.e. for 
different values of the constant magnetic field) /i = (-0.75, -0.7, -0.5, -0.2, 0.0, 0.1, 0.2, 0.4). The value of v has been fixed at 
v = 0.075. The initial values of the position vectors (measured in A = c/ui) and velocity vectors (measured in c) are (60,0,0) 
and (0.001, 0, 0.001), respectively. The size of the area shown in these plots is 100A x 100A and the time lapse is 600/cj. The 
third plot shows the exponential growth of the distance from the vortex line characteristic of the unstable regime. 



y y y y 




FIG. 5: The trajectories of a particle obtained for the same initial conditions as in Fig. 2] In these plots the value of is fixed 
at /i = —0.2 but the strength of the field amplitude varies as follows: v — (—0.1, —0.05, 0.05, 0.1). 



and magnetic field vectors read 



Bu ' 

cB c 
Buj 



x cos i? + y sin •& — 
x sin ■& — y cos $ - 

' y cos "d — x sin $ 
y sin § + x cos i? 



X\ COSt?i - 

x\ sin??i 4 


- yi cos??i 

- yi sini?i • 




- yi sintfA 

•yi cosi?i , (45a) 

+ x\ sint?!^ 

- xx cost?i I , (45b) 



where $ — u(t — z/c) and #i = u)\{t — z/c). Thus, the 
vortex line forms a helix rotating with the angular fre- 
quency lo — u}\. The projection of this helix on the xy 
plane forms a circle. In Fig.|Hlwe show the motion of the 
particle dragged by such a helical vortex line. The values 
of the dimensionless parameters [i = 0.1 and v = 0.12 
that determine the character of the motion may corre- 
spond to different field and particle configurations. For 
example, they can describe an electron moving in a con- 
stant magnetic field of 1 T and in an electromagnetic 
microwave of the frequency 2.810 11 Hz and characteristic 
intensity 3.410 10 W/cm 2 . These parameters are not cho- 
sen to be very realistic (for example, the intensity is very 



high) but they were selected to enhance characteristic 
features of our solution. 

In order to obtain a more elaborate path (a generalized 
Lissajou figure) of the vortex in the xy plane, we add two 
plane waves instead of one, i.e. 

F g {r : t) =Biv(x + iy) [(x + ^e^*-^ 
-( Xl + iyx)e- l ^ {t - zjc) - [x 2 + iy 2 )e- l ^ {t - z/c) ] . (46) 

In this case, the position of the vortex as a function of t 
and z is given by the equation 

x v (z, t) + ty v (z, t) = {xx + iy 1 )e i < w - Wl W-*/°> 

+ (X 2 + iy. 2 y("-"2)(t-z/c)_ ( 47 ) 

The numerical solutions presented here in Figs. 
should leave no doubt that an electromagnetic vortex 
may easily overcome the inertia of the particle and drag 
the cyclotron orbit along. Analytic solutions of the equa- 
tions of motion for a relativistic particle in the presence 
of electromagnetic waves with moving vortices can also 
be obtained and they will be described elsewhere. 



FIG. 6: The trajectories of a particle obtained for the same initial conditions as in Fig.|3] In these plots the value of is fixed 
at fi = 0.5 but the strength of the field amplitude varies as follows v — (—0.2, —0.1, 0.1, 0.2). 




FIG. 7: The transport of an orbit by a moving vortex line. 
On top of the average motion with a constant velocity in the 
x direction there are additional oscillations, clearly seen after 
a tenfold enlargement in the upper part of the figure. 



V. QUANTUM MECHANICS 

We shall first describe the motion of a quantum particle 
in the combination of an electromagnetic wave with a 
vortex line [Eq.JJJJ] and a constant magnetic field in a 
fully relativistic case. The nonrelativistic limit will be 
obtained by a simple reinterpretation of the parameters. 
Also, the incorporation of spin does not introduce any 
significant complications. Therefore, we shall start with 
the Dirac equation. 

A. Relativistic quantum mechanics 



FIG. 8: The transport of an orbit by a vortex line sweeping 
out the surface of a cylinder with the radius 10 A. The field 
parameters are pi = 0.1 and v = 0.12. The actual particle 
trajectory (upper plot) is compared with the path that traces 
the points of the intersection of the vortex line with the xy 
plane defined by the equation z = £[t] (lower plot). The 
period of the circular motion was chosen as 1000 times the 
period T of the carrier wave and the time lapse for the tra- 
jectory shown in this figure is equal to four full periods of the 
circular motion. The distances are measured in units of the 
wavelength A = c/ui of the carrier wave. Initially the particle 
was placed at the position of the vortex (£o = 10, rjo = 0) 
with zero velocity. A slow drift in the z direction is due to 
relativistic effects. 



The Dirac equation for a particle moving in the pres- 
ence of an electromagnetic field described by a vector 
potential A(r, t) can be written in the form 

ihdt^ = [ca-{-ihV - eA) + fimc 2 ] *. (48) 

In our case, the vector potential can be chosen as 



-Bg(r,t)-yB /2 
Air J) = j Bf(r,t)+xB /2 



(49) 



The quantum counterpart of the classical symmetry 
transformation (|13fl . namely 



-a, x 



-OJQ, * 



1 + a x a y 
V2 



*, (50) 



leaves the Dirac equation invariant. Therefore, we may 
again restrict ourselves to the case a = 1 and obtain the 
solution for a = — 1 by applying the transformation 150|l . 
It is convenient to choose the Dirac representation |24j of 
on and (3 but with the matrices a z and (3 interchanged, 
i.e. 

a x = p x a Xl a y = p x a y , a z — p z , (3 = p x a z , (51) 

where pi and Ui are two sets of Pauli matrices acting in- 
dependently. Next, we split the four component bispinor 
into two-component wave functions with the help of two 
projection operators P± = (1 ± a z )/2. In the represen- 
tation H51f) . the matrix a z is diagonal and the splitting 
amounts to taking the upper and lower components of 
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and y components undergo a rotation, for example, 
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FIG. 9: The transport of an orbit by a vortex line sweeping 
out a cylindrical surface with the cross section in the form of 
a trefoil. This figure shows only the projection of the trajec- 
tory on the xy plane. The two frequencies in the formula 14611 
are 1.0005 a; and 0.999 u. Since the two detunings are com- 
mensurate, the trajectory (right) of the vortex line in the xy 
plane is closed. The parameters characterizing the position 
of the vortex are x\ = 5, y% = 0, X2 = 0, and j/2 = —10. The 
remaining parameters are fj, — 10 and v — 1.5. The distances 
are measured in units of the wavelength X — c/cu of the car- 
rier wave. Initially the particle is placed on the vortex line at 
£o = 5 and rjo = — 10 with zero velocity. The size of the area 
in these figures is 20A x 20X. 



the Dirac bispinor Therefore, we may write 

* = P + ¥ + P_*=f* + ) + f^°_) (52) 
and the Dirac equation l|48|) takes on the form 



2iHd + % + = 



(53a) 
(53b) 



where d± = d/dt± = (d t ± cd z )/2 and 

h = c [mca z — (t± ■ (ihV + eA)] . (54) 

From now on, all capitalized Greek letters denote two- 
component spinor wave functions. In the next step, fol- 
lowing the procedure of Ref . , we perform the transfor- 
mation to the rotating frame to eliminate the dependence 
of the potential on t — z /c. To this end, we make the fol- 
lowing substitutions in Eqs. (f 531) 

(55) 



where M z is the z component of the total angular mo- 
mentum 

M x = L z + ~o z = ~(xd y - yd x ) + -a z . (56) 

Under the transformation 1)55(1 . the z components of all 
vector operators, r, V, and cr do not change while the x 



U'xU = xcos{ujt-/2) — 2/sin(arf_/2), 
U^yU = xsm(uit-/2) + ycos(wi_/2). 

Applying these transformations to both Eqs. 
(I53bl> we obtain 



{2ihd- 



2ihd+^ A 



where 



(57a) 
(57b) 

aj and 

(58a) 
(58b) 

(59) 



ho — c [mca z — er± ■ (ihV + cAq)] , 

and Aq is the vector potential 149|) evaluated at i_ = 0, 

A = ((B-B Q /2)y,(B + B /2)x,0). (60) 

Since now there is no explicit dependence on t and z, 
we may seek solutions of Eqs. (|58|) by separating these 
variables 



^ ± (x,y 1 z 1 t) = e-^ Et -^ h ^ ± (x,y), 



(61) 



where E and p z are the separation constants. The equa- 
tions for i&±(%, y) have the form 

(E- Pz c)<S>+ = hoi-, (62a) 
(E+p z c + wM z )$- = h $ + , (62b) 

and upon the elimination of $ + , we obtain the following 
eigenvalue equation for the transverse motion 



(ihW + eA ) 



n 



M z = (63) 



2m 2 2 

where the energy of the transverse motion is defined as 

E 2 - m 2 c 4 - {p z c^ 2 



E\ 



2 A 

m c - 
2mc 2 



(64) 



and the effective frequency il = £lu = (E — p z c)u> /mc 2 is 
the same as in the classical theory [cf. Eq. 1)15(1 ], 

We shall now determine the eigenvalues and the eigen- 
functions of the Hamiltonian appearing in 1(631) . Upon 
the substitution of the explicit form ((60(1 of the vector 
potential Aq, we obtain the following equation 



h 2 m / uj \ 2 2 m / uj 

2^ a± + ^{ UJc+ y) x + 2ry 



+ihuj c (xd y + yd x ) 



Wq + f 



L0 



n 



L, 



$_=£ ± $_, (65) 



where L z is the z component of the orbital angular mo- 
mentum. This equation, by a substitution 



(66) 
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is transformed to the final form 



where 



OJ c VL ) X + 



n . ml ojP. 

— A ± H — 

2m 2 V 4 



UJ{) + il f ( LOq Q 

' T ' 7 1 



m I w. 



f + ^ y 2 



(67) 



The Hamiltonian in this equation differs only by the spin 
term from the classical Hamiltonian J22J). Since the spin 
part contributes only an energy shift, we may use the 
classical amplitudes of normal modes [Ea.l|31|)] to con- 
struct the annihilation operators 



a+ 



{HN+y/2 



•'++ 



1 



+iet- 



-ist+_ I ~d x 
i 



t— -a 



rd x 



i + f-i 

f - 

1 + fi 



1 



-y 



(68a) 



(68b) 



The normalization of these operators is different from 
their classical counterparts by a factor of l/Vh to secure 
the proper commutation relations. The quantum Hamil- 
tonian, appearing in Eq. I|67|) , expressed in terms of the 
creation and annihilation operators is (without the spin 
term) 



H = /i£l + a^a + — sgn(l/2 + /i)/i£l_dla- 



sgn(l/2 + fj,) 



m- 



(69) 



where we have used the relations £±+i±- = 2|1 + /i|£!±. 
It differs from the classical expression ((35(1 by the usual 
zero-point energy terms. 

Now, we shall determine a Gaussian wave function an- 
nihilated by the operators a±. We may do so by substi- 
tuting a Gaussian wave function in the form 

i'o = cxp(-q x x 2 /2 - q y y 2 /2 - iqxy) 

into the equations 

a + ipo — 0, a_-0o = 



(70) 



(71) 

and find q x ,q y , and q. By comparing the coefficients at 
x and y in (|71|) , we obtain the following set of four linear 
equations for the three parameters q x , q y and q 

q x \l + fi\ i+_ - q(l + fi) t ++ = mQs-t ++ /H, (72a) 

q y \l + fi\t ++ +q(l+ fi)t+- =mfls + t + -/H, (72b) 

q x \l + n\ <_+ + q(l + n)t = mQs+t /h, (72c) 

q y \l + n\ t - q(l + /i) t- + = mOs_t_+/fi. (72d) 

There exists a solution of these equations and it reads 

q x = zu(s + + s-)t ++ t — , (73a) 
q y = m(s + + s_)< + _t_ + , (73b) 
q = em(s + t^ — s_t ++ t |_), (73c) 



h\l + +t+_i__)' 



(74) 



The parameters q x ,q y , and q can also be expressed as the 
following explicit functions of /i and v 



qx 



1 + 2/i - Av 
mfl 



(75a) 



% = 



x y/(l + 2/i - 4i/) [(1 + 2/i)(l + v) 2 - 8^ 2 - u], (75b) 



8hv 2 



(1 + 2/i-u), 



(75c) 



where u = sgn(l/2 + fi)y/(l + 2/i) 2 - 16^ 2 . 

The Gaussian wave function l|70() describes an analog of 
the ground state of the system. This will not always be a 
true ground state because for > —1/2 the Hamiltonian 
is not positive definite. Having found the state tpo we 
may generate a complete set ip mn of "excited" states by 
acting on tJjq with powers of the creation operators 



:na — a + a + ^o- 



(76) 



The wave functions representing these states are Gaus- 
sians [Eq. (|70(l ] multiplied by the polynomials in the vari- 
ables x and y of the (m -I- n)th order. All wave func- 
tions [Eg . 176( 1] are localized around the vortex line. A 
complete set of two-component wave functions can be 
obtained from the one-component functions ip m n by at- 
taching two independent two-component spinors u± 



ip mn u±. 



(77) 



Choosing the spinors u± in the form u + = (1,0) and 
u_ = (0, 1) we obtain for each choice of m and n two so- 
lutions of the eigenvalue equation l|67l) of the same shape 
and differing only in the energy eigenvalues. 

Having reduced the problem to a two-dimensional har- 
monic oscillator, we may use the whole arsenal of tools 
available in this case. We may combine the classical 
trajectories with the solutions of the wave equations 
and construct quantum counterparts of classical motions. 
Namely, with each pair made of a classical solution of 
the equations of motion described by the Hamiltonian 
[Eg. 123(1 ] and a localized wave function [Eq. 1(76(1 ] we 
may associate a new solution of the wave equation in the 
form of a displaced wave function [2^, [2(| [27j represent- 
ing a localized wave packet moving along the classical 
trajectory. These states correspond to coherent states of 
quantum optics. We may also introduce the analogs of 
optical squeezed states described by Gaussian wave func- 
tions whose shapes are changing during the time evolu- 
tion. In the next subsection we shall use the Gaussian 
states to prove that the transport of the orbits by moving 
vortices, clearly seen in the classical case, carries over to 
the quantum case. 
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B. Nonrelativistic quantum mechanics 

The eigenvalue equation l|63|) does not contain the 
speed of light. Therefore, one may suspect that it co- 
incides with the energy eigenvalue problem for a nonrel- 
ativistic charged particle in a coordinate frame rotating 
with the angular frequency Q/2. We may confirm this ob- 
servation by neglecting the relativistic corrections from 
the very beginning. This amounts to using the nonrela- 
tivistic Schrodinger-Pauli equation instead of the Dirac 
equation 



'ih\7 + eA nr ) 2 eh 



2m 



2m 



<7 Bn * 



(78) 



where in the spirit of the nonrelativistic approximation 
we have neglected retardation in the expression (|49|l for 
the vector potential 

B(y cosuit — x sin tot) — yBo/2 
A nr (r,t)= ( B(x cos u;t + y smut) + xB /2 ] (7!)) 

Note that in the Schrodinger-Pauli equation the mag- 
netic moment is coupled only to Bo because the magnetic 
field generated by the vector potential (|79|l has only the 
constant part. After the transformation to the rotating 



frame we obtain the following eigenvalue equation 
'(ih\J + eA Q ) 2 huj 



2m 



., a g --M s )9 = m. (80) 



It differs from its relativistic counterpart [Eq. (|63[) ] only 
by the replacement of the relativistic parameters f2 and 
E± by their nonrelativistic limits oj and E. Thus, the 
mathematical solutions of the eigenvalue problem ob- 
tained in the full relativistic theory can be used, without 
any modifications, in the nonrelativistic case. This does 
not mean that the relativistic corrections vanish. They 
are all contained in a difference between the relativistic 
parameters (f2, E±) and the nonrelativistic parameters 
(uj,E) characterizing these solutions. 

We shall now apply the time dependent Schrodinger- 
Pauli equation to the case of a moving vortex to show 
that the motion of a quantum wave packet corresponds 
to its classical counterparts. It is convenient to work 
in the gauge in which the electromagnetic field of the 
wave is described by the scalar potential. This cannot 
be achieved for the full electromagnetic field [Eq. 
but it can be done in the nonrelativistic approximation, 
when the retardation effects and the magnetic field of the 
wave are neglected. In the case of a vortex moving on a 
circle, the electric field [Eq. (|45a|l ] is obtained from the 
following scalar potential 



V(r,t) = -eBuj 



2 2 

x- — y z 



cos(uit) + xysin(ujt) — x{x\ cos(wii) + y\ sin(u;ii)) — y(—yi cos(wii) + X\ sin(ojii)) I . (81) 



Since the variable z does not appear in the potential, we can separate out the z-dependent part of the wave function 
and consider only the following Schrodinger-Pauli equation describing the motion in the xy plane 



ihd t ^{x,y,t) 



■5-A ± + — -^(x 2 + y z ) - -hw (xdy - yd x ) + V{x, y, t) 
2m 8 2i 



Hujq 



(82) 



Upon the transformation to a frame rotating with the angular velocity u>/2, the time dependence in the quadratic 
part of the Hamiltonian disappears and we obtain an equation for a stationary driven oscillator 



ihd t ^f(x,y,t) = 



h 2 A mu> 2 9 u>q + uj ~ x 2 — y 2 

+ — H* 2 + V 2 ) - ^7^L Z + u:^—^ 
2m 8 2 2 



. ' Hun hul \ 



*(*,!/,*), (83) 



where E\ (t) is the part of the electric field in l|45a() with 
the frequency lo\. In the absence of the driving force 
E\ (t) , the Gaussian solution of this equation is given by 
{7UJ. In order to construct the solution of this equa- 
tion that corresponds to the classical trajectory shown in 
Fig. |H1 we shall employ a procedure that enables one to 
obtain from any given solution ipo(r, t) of the Schrodinger 
equation with a quadratic Hamiltonian a family of alter- 
nate solutions. These solutions are labeled by all possible 
classical trajectories and are obtained by "space-shifting" 
and "phase-shifting" the original solution according to 



the formula 

^ } (r, t) = e-^/^^e-^^Mr, t), (84) 

where £*(<) and iTi(t) are the classical positions and mo- 
menta and the phase S(t) is the integrated classical ac- 
tion. This construction is based on a generalization of the 
Dobson's harmonic potential theorem |2(|. In the Ap- 
pendix we give, for completeness, a simple proof of this 
theorem for the most general quadratic Hamiltonian. We 
shall consider the simplest case, when ip in the rotating 
frame is given by the Gaussian wave function 1)70(1. In 
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FIG. 10: The transport of a quantum mechanical wave packet in the xy plane by a moving vortex line. Four snapshots of 
the probability distribution show the wave packet moving along the classical trajectory (superimposed for comparison). The 
parameters of the classical trajectory are the same as in Fig. [HI The contour lines enclose the areas where the particle is found 
with the probability 0.1, 0.5, 0.9, and 0.99, correspondingly. Note that due to the influence of the vortex wave the wave packet 
is squeezed and rotating. 



the laboratory frame, the corresponding wave function is 
obtained by the inverse transformation of the variables 
as compared to Ea. (|57|l and it has the form 

ip L (r,t) — e -9^[ a:cos (" t / 2 )+a sin ( wt / 2 )] 2 / 2 



xe 



-q y [y cos(aji/2) — x sin(ujt/2)] 2 /2 



xe 



-iq[x cos(u;i/2)+y s\n{ujt / 2)][y cos{ujt/2)—x sm(oji/2)] (85) 

The new solution corresponding to the classical trajec- 
tory depicted in Fig. [S] obtained according to the pre- 
scription (|84|l is 

7r )(r,i) = e - i S(t)/h e ~q m [xcos(u J t/2)+ysm(u J t/2)] 2 /2 
y^ e ~1y[v cos(wt/2)— £csin(wt/2)] 2 /2 
^ g — iq[x cos(uJi/2)+y sin(ujt/2)] [y cos(uji/2) — x sin(uJi/2)] (86) 

In Fig. llUl we show the time evolution of this wave packet, 
when the center of mass motion is described by the non- 
relativistic limit of the same classical trajectory that is 
shown in Fig. [S] The Landau orbit is transported along 
a circle by the moving vortex. This wave function de- 
scribes an electron localized with probability 0.99 inside 
the largest ellipse, transported in the constant magnetic 
field of by an electromagnetic wave with a moving vortex. 



netic field can be controlled by electromagnetic waves 
with embedded vortex lines. These orbits arc pinned 
down in the vicinity of the vortex line and are dragged 
along when the vortex line moves in the plane perpen- 
dicular to the magnetic field. We analyzed this behavior 
with the use of our new analytic solutions of the Lorentz 
equations of motion in the classical case and the Dirac 
(or Schrodinger) wave equation in the quantum case. We 
show that the trapping of orbits by the vortex in the 
classical case has its counterpart in the form of Gaussian 
shape of the wave functions localized in the vicinity of the 
vortex. The effective reduction of the dynamics of our 
system to that of a two-dimensional harmonic oscillator 
makes it possible to give a complete solution using many 
tools developed in the past for oscillators. In particular, 
with the use of creation operators we have constructed 
a complete set of states. Coherent states, the closest 
analogs of classical orbits, are also easily constructed. 
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APPENDIX A: SOLUTIONS FOR A GENERAL 
QUADRATIC HAMILTONIAN 

In this appendix we show how to construct the new 
solutions of the Schrodinger equation that were used to 
describe the transport of Landau orbits. Let us consider 
the case of the most general quadratic Hamiltonian 



1 



H(x\ Pi ,t) = -a 11 {t)it^3 + 



1 

fi{t)C 



(Al) 



where the summation over repeated indices is under- 
stood. The Hamiltonian equations of motion have the 
form 

£ k (t) = a k > (t)n 3 + efit)? (*) + ff*(t), (A2a) 
= -b kj (t)e - c£{t)*j(t) + f k (t). (A2b) 

The Schrodinger equation for this system is (ft = 1) 



id t ip(x l , t) 



-a? j (t)d i d j + -b ij (t)x i xi 



- icP (t)x l d J: -Mt)x l - igitfdi + -c/Jv^V)- (A3) 

The last term on the right hand side is added to make 
the quantum Hamiltonian a Hcrmitian operator. For this 
system, the following statement holds. For every solu- 
tion ijjo(x l ,t) of the Schrodinger equation (|A3I) without 
external forces {fi — — g 1 ) there corresponds a family 
of solutions ip£ n (x l , t) of the complete equation i|A3|l ob- 
tained by space-shifting and phase-shifting the original 



solution, namely, 

^(x\ t) = e tS ^e^ x 'M^ ~ ?(t), t) (A4) 

where (£ l (t), Ti(i)) is any solution of the classical equa- 
tions of motion \A2\ , To prove this statement we sub- 
stitute into i|A3(l the new wave function written in the 
form 

iP in {x\t) = e- iS ^U v U £ M^,t), 

U T = e l ^ (t)x \ Us = e-?® 8 * (A5) 

and we multiply the whole equation from the left by U~ , 
U^ , and e~ tS ^\ Next, we use the relations 

U-'-dkU* = d k + in k (t), U^x k U^ = x k + £ k (t), 

e iS(t)g te -iS(t) = g t _ idtS (t), 



U- l d t U v = d t + iMt)x\ U^dtUt = 8t- €(t)di (A6) 



and we collect all the terms linear in x % and di and in- 
dependent of x % and di. The linear terms cancel due to 
the classical equations of motion i|A2(l and the indepen- 
dent terms cancel if the phase S(t) is made equal to the 
following classical action integral 

s(t) = J dt - \h,ee + g 1 ^ . (A7) 

All remaining terms reduce to the Schrodinger equation 
satisfied by the wave function tpoi^ 1 , t). 
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